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We study the ground state of a clean two-dimensional electron liquid in a weak magnetic field 
where N 1 lower Landau levels are completely filled and the upper level is partially filled. It is 
shown that the electrons at the upper Landau level form domains with filling factor equal to one 
and zero. The domains alternate with a spatial period of order of the cyclotron radius, which is 
much larger than the interparticle distance at the upper Landau level. The one-particle density of 
states, which can be probed by tunneling experiments, is shown to have a gap linearly dependent 
on the magnetic field in the limit of large N. 
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The nature of the ground state of an interacting two- 
dimensional (2D) electron gas in magnetic field has at- 
tracted much attention. The studies have been focused 
mostly on the case of very strong magnetic fields where 
only the lowest Landau level (LL) is occupied, so that 
the filling factor v — k\l 2 does not exceed unity (here kp 
is the Fermi wave- vector of the 2D gas in zero magnetic 
field and I is the magnetic length, I 2 = H/muj c ). The 
physics at the lowest LL turned out to be so rich that, 
perhaps, only at v — 1 the ground state has a simple 
structure. Namely, it corresponds to one fully occupied 
spin subband of the lowest LL. The charge density in such 
a state is uniform. The case of a partial filling, v < 1, 
is much more interesting. Using the Hartree-Fock (HF) 
approximation, Fukuyama et al. |]J found that a uniform 
uncorrelated spin-polarized electron liquid (UEL) is un- 
stable against the formation of a charge density wave 
(CDW) at wave- vectors larger than 0.79 The opti- 
mal CDW period was later found to coincide with that 
of the classical Wigner crystal (WC) 0. 

Subsequently, however, it turned out that non-HF trial 
states suggested by Laughlin || f° r v = 1/3, 1/5 to ex- 
plain the fractional quantum Hall effect are by a few per- 
cent lower in energy. The Laughlin states were further 
interpreted in terms of an integer number of fully occu- 
pied LLs of new quasiparticles, composite fermions [Q. 
This concept was then applied to even denominator frac- 
tions Hj. Thus, although the HF approximation gives a 
rather accurate estimate of the energy, it fails to describe 
important correlations at a partially filled lowest LL. 

Recently, the requirement of the complete spin polar- 
ization in the ground state was also reconsidered, ft 
was found that a partially filled lowest LL may contain 
skyrmions ||. 

In this Letter we consider the case of weak magnetic 
fields or high LL numbers N. There is growing evidence 
from analytical and numerical calculations that fractional 
states, composite fermions and skyrmions are restricted 
to the lowest and the first excited LLs (N = 0,1) only 
(see Refs. 0-||]). We will present an additional argument 



in favor of this conclusion. This point of view is also con- 
sistent with the experiment because none of those struc- 
tures has been observed for N > 1. 

Before we proceed to the main subject of the paper, a 
partially filled upper LL, note that we can use the con- 
cept of LLs only if the electron-electron interactions do 
not destroy the Landau quantization. For weak magnetic 
fields where the cyclotron gap hui c is small, this is far 
from being evident. To see that the LL mixing is indeed 
small one has to calculate the interaction energy per par- 
ticle at the upper LL and verify that its absolute value is 
much smaller than Huj c . The largest value of the interac- 
tion energy is attained at v = 2N + 1 where the electron 
density at the upper LL is the largest. The interaction 
energy per particle is equal to — ^E cx , where E cx is the 
exchange-enhanced gap for the spin-flip excitations |ic| l 
at v = 2N + 1 (it determines, e.g., the activation energy 
between spin-resolved quantum Hall resistivity peaks). 
Aleiner and Glazman (AG) || calculated E cx to be 
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where Eh_ is the "hydrodynamic" term (see Ref. | fTl|| ) 
given by |1 



Eh = huj, 
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The parameter r s entering these formulae is defined by 
r s = V2/kpaB, ob = h 2 n/me 2 being the effective Bohr 
radius. In realistic samples r s ~ 1 but even at such 
r s the ratio E ex /Hu> c is still rather small. Therefore, 
even at weak magnetic fields the cyclotron motion is pre- 
served and the mixing of the LLs is small. Note that the 
first term in E ex linearly depends on the magnetic field 
whereas Eh_ has an approximately quadratic dependence. 

Since we chose to rely on the HF approximation, a 
natural turn of thought is to consider a WC-type state 
whose wave-function is given by P]l3[] 
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where |Ojv) stands for N completely filled LLs and c^ R 
is the creation operator for a certain one-particle state, 
called a coherent state [|l4| . The modulus of the coherent 
state wave-function is not small only within a distance I 
off the classical cyclotron orbit with the center at the 
point R and radius R c = kpl 2 . In the HF WC state Ri 
coincide with the sites of a triangular lattice with density 
vn / (2tt1 2 ) , where vn = v—2N. From now on we consider 
only j//v < \, which suffices because of the electron- hole 
symmetry. 

When vn is small, vn <C 1/-ZV, the cyclotron orbits at 
neighboring lattice sites do not overlap, and the concept 
of the WC is natural. However, this concept was ap- 
plied for overlapping orbits as well. According to AG, at 
N > r\T 2 > 1 and not too small v N , v N » l/(Nr 2 s ), the 
cohesive energy of the WC, i.e., the energy per particle 
at the upper LL with respect to that in the UEL of the 
same density, is given by pa 
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Assuming that the WC is the ground state, AG found 
that the one-particle density of states (DOS) consists of 
two narrow peaks separated by the pseudogap E s = E\ x 
(see also Ref. 111]]). In the limit of large N, both E g and 
l-^coh I are mucn smaller than E ex , and so AG concluded 
that there are two different scales for spin and charge 
excitations. 

In this Letter we claim that for i/jy l/(Nr 2 ) the 
ground state is not the WC, but another type of a CDW 
whose period is of order R c . In contrast to the lowest 
LL, the optimal CDW period is much larger than the 
average distance between the electrons at the upper LL. 
The cohesive energy of the CDW has the scale E cx and 
is given by 
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where /(vn) ~ 0.03 at vn = h and /(fjv) oc vn at 
l/(Nr 2 ) < v N < §. The DOS consists of two peaks 
(van Hove singularities) at the edges of the spectrum, 
the distance between them for vn ~ | being equal to 
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Hence, we claim that all the important properties of N-th 
LL are determined by the single scale, E ex . 

Let us compare E™£ and E^ w . The "hydrody- 
namic" term is the same in both. Hence, one has to com- 
pare only the remaining terms. It is easy to see that the 
CDW state wins over the WC provided v N > l/(Nr 2 ). 

Our CDW state can be roughly approximated by a 
state (||), with Ri forming patterns shown in Fig. |l|. The 
aggregation of many particles in large domains of size 



R c allows the system to achieve a lower value of the ex- 
change energy. At the same time due to the fact that 
the domain separation is chosen according to the spe- 
cial ring-like shape of the wave-functions at the upper 
LL, the actual charge density variations are not too large 
(of order 20%). Hence, the increase in the Hartree en- 
ergy due to the domain formation is small. According 
to our numerical simulations for N — 5 and r s = 0.5, 
at vn > 0.3 the optimal CDW has a "stripe" structure 
(Fig. 0a). At v N < 0.3 a "bubble" pattern (Fig. |l]b) 
wins. The distance between the "bubbles" in this pat- 
tern is of order R c and remains approximately the same 
as vn decreases. Correspondingly, their diameter is given 
by ~ R C y/i7N. At vn ~ 1/N where it becomes of order I, 
the "bubbles" consist of single electrons, i.e., the CDW 
state becomes indistinguishable from the WC (Fig. |l|c). 
With further decrease in vn, the distance between the 
electrons increases. 
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FIG. 1. CDW patterns, (a) Stripe pattern, (b) Bubble 
pattern, (c) WC. One cyclotron orbit is shown. 

At this point we would like to address the issue of 
the fractional states at high LLs. We believe that at 
vn S> 1/N, the fractional states can not compete with 
the CDW state. Indeed, the CDW state has a very 
low energy because of the correlations in the positions 
of the guiding centers on the length scale R c , which is 
the largest length scale in a not too dilute system. In the 
fractional states, just like in the WC, these correlations 
have the length scale I. Based on the example of the 
WC, it seems very plausible that the correlations of this 
type are much less effective. On the other hand, there 
is no doubt that at vn <C l/N the WC is the ground 
state. This leaves only a narrow window in the vicinity 
of vn = l/N, where the fractional states may or may not 
appear. 

The novel ground state enables us to explain two in- 
teresting experimental findings. One is the magnitude of 
a pseudogap in the tunneling DOS, first observed in ex- 
periments on single quantum well JH| and, recently, on 
double quantum well high- mobility GaAs systems |l7],|l8| . 
The pseudogap -Etun appears to be linear in magnetic 
field for 1 < N < 4 |l8). Theoretically, the pseudogap is 
given by E tun — 2E g . The additional factor of two arises 
because the tunneling DOS is the convolution of the DOS 
of the two wells. For the parameters of Ref. |lq] Eq. (|^) 
leads to -Etun ~ 0.52fto; c , which compares favorably with 
the experimental value of 0A5huj c |T§| ]. In the experimen- 
tal range of parameters the "hydrodynamic" term dom- 
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inates, and our result is only by 35% larger than that 
of AG, 2E\ 1 . However, in the limit JV > 1 we predict 
a much wider pseudogap with a linear instead of an ap- 
proximately quadratic dependence on the magnetic field. 
Note that even for 1 < N < 4 the dependence, which we 
predict, is not much different from the linear one. Re- 
cently, Levitov and Shytov |L9| obtained an expression 
for Stun similar but not identical to ours without study- 
ing the ground state of the system. We believe that only 
the CDW ground state can justify this type of expression. 

Another important application of the proposed picture 
concerns with the conductivity peak width of the integer 
quantum Hall effect in high- mobility structures where the 
disorder is believed to be long-range. A semiclassical elec- 
trostatic model of Efros Q predicts that the electron liq- 
uid is compressible in a large fraction of the sample area. 
If compressible liquid is considered to be metallic, then 
the conductivity peaks are necessarily wide ]|(|, which 
is indeed observed at relatively high temperatures |pl| . 
However, it is well-known that at low temperatures the 
peaks are narrow (see, e.g., Ref. ^2j), which may result 
from the pinning of the compressible liquid p3] | . The fine 
CDW structure of the compressible liquid (Fig. makes 
such a pinning possible even though the disorder is long- 
range. Note that although the pinning prohibits sliding 
of the CDW as a whole, the current can still flow along 
the boundaries of the filled and empty regions (the "bulk 
edge states"). Precisely at z/jv = |, the bulk edge states 
form a percolating network, which leads to a narrow peak 
in conductivity with, in certain models p3 ], a universal 
height 0.5e 2 /h. 

We start our analysis by writing down the HF cohe- 
sive energy of the electrons at the upper partially filled 
LL (cf. Refs. @@), 

" L ^uMq)\Mq)\ 2 - 



TP 

E. 



CDW 
coh 



2VN +n 
q^O 



(7) 



Here and below we use tilde for Fourier transformed 
quantities, L is the size of the system, = (2tt1 2 )^ 1 , 
and A(r) is the CDW order parameter. It is propor- 
tional to the guiding center density at the point r. For 
instance, the WC corresponds to A(r) in the form 0] 
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The HF interaction potential mhf(<z) entering Eq. (0) 
is the difference of the direct and the exchange terms, 
unF(q) — UB_(q) — u ex (q), which are further defined by 
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Ljy being the Laguerre polynomial. Following Ref. || 
(see also Ref. pj|), the screening by the lower LLs is ex- 



plicitly taken into account with the help of the dielectric 
constant 
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From Eqs. (P,^o|) an asymptotic expression for uhf(<?) 
can be derived, 
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We want to find the distribution of the guiding center 
density A(x,y) that minimizes the energy. Generally, 
this is a non-trivial problem because the HF equations 
have to be solved self-consistently. However, if the CDW 
is unidirectional, i.e., if A(x,y) does not depend on y : 
the self-consistency condition is simply 



A(x) - e[-eu F (x)]/L 2 , 
£hf(x) = ^2n L u HF (q)A(qx)e lq 



(13) 
(14) 



where £hf(^) is the HF self-energy, 0(x) is the step func- 
tion, and a: is a unit vector in the x direction. The mean- 
ing of this condition is that the states above the Fermi- 
level are empty and below the Fermi-level arc filled. 

For N > the Hartree potential un(q) necessarily has 
zeros due to the factor F(q) containing the Laguerre poly- 
nomial [Eqs. (PJiOl)]. The first zero, qo, is approximately 
given by qo ~ 2A/R C . The exchange potential is always 
positive; hence, there exist q where the total HF poten- 
tial Uhf is negative [in agreement with Eq. (|l2|)]. This 
leads to the CDW instability because the energy can be 
reduced by creating a perturbation at any of such wave- 
vectors (cf. Ref. @). 

In the parameter range 0.06 < r s < 1 and N < 50 well 
covering all cases of practical interest, the HF potential 
is negative at all q > q and reaches its lowest value near 
q = qo (see Fig. ||). One can guess then that the lowest 
energy CDW is the one with the largest possible [under 
the conditions (|l|Jl|)] value of \A(q x)\. The CDW hav- 
ing this property consists of alternating strips A(.x) = 
and A(x) = 1/L 2 (Fig. |l|a), and non-zero A(q) are given 
by 
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provided q is an integer multiple of <7o- Our numerical 
simulations showed that at close to | this is indeed 
the correct type of the solution in the specified above 
range of r s and N, but qo should be replaced by a slightly 
smaller value of 2.3/ R c corresponding to the spatial pe- 
riod of 2.7R C . 
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FIG. 2. The Hartree, exchange, and HF potentials in 
g-space for N = 5 and r a = 0.5. 

Having established the functional form of A(x), let 
us estimate the cohesive energy B™. Performing the 
summation in Eq. (^) with the help of Eqs. ( p^| . [L5| ). 
one recovers Eq. (|J). As for the DOS, it is given by 
(^Lft)/ 71 ") |deHF/dx| . It can be verified that €kf(x) 
reaches its lowest and largest values at x = and 
x = 7r/yo, respectively. These extrema result into the 
van Hove singularities at the edges of the spectrum sepa- 
rated by the pseudogap E g — 2|e(0)|. Eq. (^) now follows 
from Eqs. (jl2]@Jl|). 

So far we discussed the unidirectional CDW, which can 
be analyzed at least partially analytically. 2D CDW pat- 
terns were studied numerically. We restricted the choice 
of A(r) to the form (||) suggested by the WC state. Re- 
call that in the WC state Ri coincide with the sites of a 
triangular lattice with density v^n^. In the simulations 
we used a different set of Ri, corresponding to the trian- 
gular lattice with the density n^. The fraction vjy of the 
total of 50 x 50 lattice sites was initially randomly pop- 
ulated and then the energy was numerically minimized 
with respect to different rearrangements of the popu- 
lated sites. The expression for the energy follows from 
Eqs. (0J): 
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gKF(Ri- Rj)(ni-Vff)(rij -un), (16) 



where 3hf(9) = ex p(~ ^Q 2 l 2 ) whf(?) & n d is the occu- 
pancy of the i-th site. In this notation the energy has a 
transparent interpretation of pairwise interaction among 
the single-electron states |cLj. In the actual simulations 
we used a slightly more accurate expression with ffHF (0 
replaced by gm(r)/[l - exp(-r 2 /2? 2 )] (cf. Refs. 
The patterns obtained from the simulations are schemat- 
ically shown in Fig. ^ and were discussed above. 

In conclusion, we have argued that the ground state 
of a partially filled upper LL in a weak magnetic field is 
a CDW with a large period of order R c . Based on this, 
we were able to explain several important experimental 
results. 
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